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ABSTRACT 

The Andromeda Galaxy (M3 1 ) is an important test case for a number of microlensing surveys 
looking for massive compact halo objects (Machos). A long-standing theoretical prediction is 
that the high inclination of the M3 1 disk should induce an asymmetry in the spatial distribu- 
tion of M3 1 Macho events, whilst the distribution of variable stars and microlensing events in 
the M3 1 disk should be symmetric. We examine the role of stars in the M3 1 visible spheroid 
as both lenses and sources to microlensing events. We compute microlensing event number 
density maps and estimate pixel-lensing rates and event durations for three-component mod- 
els of M3 1 which are consistent with the observed rotation curve, surface brightness profile 
and dynamical mass estimates. Three extreme models are considered: a massive spheroid 
model; a massive disk model; and a massive Macho halo model. An important consequence 
of the spheroid is that, even if Machos are absent in M31, an asymmetric spatial signature is 
still generally expected from stellar lensing alone. The lensing of disk sources by spheroid 
stars produces an asymmetry of the same sign as that of Machos, whilst lensing by disk stars 
against spheroid sources produces an asymmetry of opposite sign. The relative mass-to-light 
ratio of the spheroid and disk populations controls which of these signatures dominates the 
overall stellar spatial distribution. We find that the inclusion of the spheroid weakens the M3 1 
Macho spatial asymmetry by about 20 — 30% over a disk-only asymmetry for the models 
considered. We also find for our models that Machos dominate over most of the far disk pro- 
vided they contribute at least ^ 25% of the halo dark matter density. This is a conservative 
limit since many stellar events are too short to be detected by present surveys. The presence 
of the spheroid also has beneficial consequences for M31 lensing surveys. The stellar spatial 
asymmetry is likely to be important in distinguishing between a spheroidal Macho halo or a 
highly flattened halo or dark matter dominated disk, since spatial asymmetries of opposing 
signs are expected in these cases. 
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1 INTRODUCTION 

In recent years the Andromeda Galaxy (M3 1) has become a key tar- 
get for microlensing experiments which are trying to de tect massive 
compact halo objects, or M achos ^uriere et al. 200ll: lcrotts et alj 
l200lHRiffeser et al. 2001). There are two reasons for this. The first 
is the development of new techniques which allow microlensing 
to be detected reliably against unresolved stellar fields, enabling 
microlensing searches to be directed toward more distant galax- 
ies. The second reason is that the evidence gathered by the first 
generation of microlensing experiments towards the Magellanic 
Clouds remains inconclusive regardi ng the existence of Machos 
jAlcock et alfcOOOHLasserre et alJ20 00: Afonso et al. 2003). This 
is largely because the contribution of the Magellanic Clouds them- 
selves to the observed microlensing rate is highly uncertain. 

The favourably high inclination of the M31 disk (i — 77°) 
presents a very promising diagnostic for the current experiments to 



establish the presence or absence of Machos. "Crotts' noted 
that if M3 1 Machos have a spheroidal distribution the inclined disk 
should induce a spatial asymmetry in the distribution of M3 1 Ma- 
cho microlensing events, with more events being seen towards the 
far disk which lies behind a larger fraction of the halo column 
density. By contrast, foreground Milky Way Machos, stellar mi- 
crolensing events in the M3 1 disk, and variable stars mistaken for 
microlensing, are expected to be symmetrically distributed. The 
near-far asymmetry prediction therefore provides, in principle, a 
clean and simple way to distinguish M3 1 Machos from other pop- 
ulations. Several theoretical studies have exploited this to establish 
how Macho or galactic model parameters cou ld be meas ured by 
the M31 microlensing surveys jcwkTcrottifcoOO; Keri ns et alj 
l200ltlBaltz. Gvuk & Crottsl2003HKerins et all2003i) . 

The visible spheroid of M31 has the potential to com- 
plicate this simple picture. Despite superficial similarities be- 
tween our own Galaxy and M31, their visible spheroid pop- 
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ulations appear to be different in two key respects. Firstly, 
several studies indicate tfiat, as well as comprising metal-poor 
stars, the M31 spheroid contains a dominant metal-rich popula- 
tion I Holla nd, Fahlman & R icher 1996; Durrell, Harris & Pritchet 
[200 1; B ellazzini etalJl2003l) . in contrast to the exclusively metal- 
poor Galactic spheroid. Secondly, the M31 spheroid has a much 
higher luminosity density than that of our ow n Galaxy in relative 
terms CReitzel, Guhathakurta & Gould 1998). These two features 
may have a common link. Ferguson e t al. 12002) have reported the 
existence of significant stellar substructure in the outer regions of 
M31. Substructure is evident both as coherent stellar density en- 
hancements and as metallicity variations. The alignment of a stel- 
lar stream to the location of the satellite galaxy M32 indicates 
that this galaxy may be the origin of some of this substructure. 
iFerguson et alJ i2002h suggest that M32 may have once been a 
much larger galaxy and that its material has been stripped to pollute 
and beef-up the M3 1 spheroid. 

Whatever is the origin of the M31 spheroid, i ts impact on 
the exp ected microlensing signature must be assessed. iKerins et alj 
y003|) have already noted that the low surface brightness stellar 
streams are unlikely to have a significant impact for microlens- 
ing studies. The same cannot be assumed for the M31 spheroid. Its 
shape and density means that it has the potential to induce an asym- 
metry similar to that for M31 Machos, making the task of identi- 
fying a Macho population more difficult. Its stars may also act as 
sources to microlensing events arising from other lens populations, 
modifying their spatial distribution. 

In this study we investigate the likely effect of the spheroid 
population on the spatial microlensing signature towards M3 1 . Sec- 
tion |2| provides analytical arguments for asymmetry from stellar 
lenses. In section |3| we present a set of extreme but simple three- 
component galactic models for M31 which are likely to bracket 
realistic models. We compute microlensing event number density 
maps and discuss microlensing rates for these models in Section|4] 
and discuss the effect of the spheroid on Macho searches in Sec- 
tion|5l 



2 ASYMMETRY FROM STELLAR LENSES 

The predicted near-far asymmetry due to M3 1 Machos is expected 
to occur because of the larger halo column density lying i n front 
of the far disk compared to that in front of the near disk iCrottsI 
Il992h . We may therefore anticipate a similar signature due to M31 
spheroid stars acting as lenses against disk sources. However, in 
this case there is also a reverse asymmetry due to lenses in the disk 
microlensing background spheroid stars. Which of these two asym- 
metry signatures might we expect to dominate the spatial signature 
of stellar lensing events? We can answer this question by consider- 
ing the number density of stellar lensing events in the two compo- 
nents. The event number density along a given line of sight scales 
with lens (source) distance Di (Ds) as 



^vpi{Di)dDips{D,)(^^^ 



dDs, (1) 



where p\ (ps) is the lens (source) mass density, {M/ I/)s is the mass- 
to-light ratio of the source population and v is the mean pairwise 
velocity of lenses and sources. Equation is proportional to the 
rate of microlensing for sources at Ds, multiplied by the source 
number density. 

For a fixed source luminosity density the number of sources 
depends upon the typical source luminosity Ls. However, the ex- 



pected number of pixel-lensing events turns out to be independent 
of source luminosity to first order. The faintest detectable pixel- 
lensing events require a minimum source magnification Amin oc 
LjT^ in order to be seen against the local surface brightness flux. In 
pixel lensing the large bulk of observed events involve highly mag- 
nified sources, for which Amin — Wmax ™d ""max is the largest 
impact parameter between lens and source which still provides a 
detectable event. Since the pixel-lensing rate per source Fp oc Umax 
it follows that Fp oc Ls- So, even though the number of potential 
sources A's oc L^^ for a fixed surface brightness, the total number 
of pixel-lensing events A'^ — N^Tp does not depend on Ls to first 
order. 

For lines of sight passing through the major axis of M31 the 
disk stars are, to a good approximation, at a fixed distance Dm-h 
from the observer. So, in the case of spheroid lenses and disk 
sources we have psdDs = Edisk5(Ds — Dm3i), where Edisk is the 
disk column density and 5{x) is the Dirac delta function, pi = psph 
is the lens density distribution and [M / L)b = (M/I/)disk is 
the source star mass-to-light ratio. The event number density for 
spheroid lenses and disk sources for lines of sight passing through 
the M3 1 major axis is therefore 
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In the case of disk lenses and spheroid sources we have ps = Pspii, 
pidA = Edisk<5(Ai3i - A) and (M/L)s = (M/L)sph, imply- 
ing an event number density 



diVdisk oc E, 



iiskf-V' 

\ L J spii 



Dm3i(A- A 



-Psph(A)dA.(3) 



spli V Ds 

In both equation J2j and J3} above we have dropped the factor v 
since the symmetry between the two lens populations we are con- 
sidering permits us to treat it as a constant. Given the spatial sym- 
metry of the disk and spheroid populations for lines of sight pass- 
ing through the M3I major axis we can compare directly the rate 
contribution dA'gph for spheroid lenses at Di — Dm3i — z to the 
contribution dA'disk for sources at A = Dm3i + z: 

dNsph _ 1 (A//L)sph 



dA^di, 



-Dm3i (M/L)disk 
(M/L)sph 



\/ {Dm-m — z){Dm31 + z) 



{z < Dm3i) 



(4) 



(A//L)disk 

This simple formula tells us that the direction of asymmetry for 
stellar lenses is controlled by the mass-to-light ratio of the disk and 
spheroid populations. If M / L for the spheroid exceeds that for the 
disk then one expects spheroid lenses to dominate number counts 
on the major axis, and so the overall stellar asymmetry should be in 
the same direction as that of M31 Machos, i.e. more events towards 
the far disk. If the disk has a larger M/L then disk lenses should 
dominate and we can expect the stellar asymmetry to run counter 
to that of M3 1 Machos. 



3 EXTREME MODELS OF ANDROMEDA 

In assessing the likely contribution of the visible spheroid popula- 
tion we choose to set up three extreme models which are likely to 
bracket realistic models. We are not trying to present the models as 
likely descriptions of the M31 structure, rather we wish to ensure 
that the real signature lies somewhere within the range of results we 
find for our extreme models. Each of our models comprises three 
components: a visible spheroid, a disk population, and a dark halo. 
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Table 1. Parameters for the three extreme M31 models. From the left the columns give, respectively, the model name, component population, 
mass of the component, its central density po . and the cut-off radius R. Additional columns give, where appropriate, the core radius a, the B- 
band mass-to-light ratio M/Lg in solar units, the disk scale length h and height H, and the flattening parameter q. Parameters for a Reference 
model, representing a less extreme "mainstream" model are also given. 



Model 


Component 


Mass (xlO^'^ Mq) 


po(M0pc-3) i?(kpc) 


a(kpc) 


M/Lg 


/i(kpc) 


H(kpc) 


q 


Massive halo 


halo 


191 


0.25 155 


2 








1 




spheroid 


4.4 


4.5 40 


1 


9 






0.6 




disk 


3.2 


0.24 40 




4.5 


6 


0.3 




Massive spheroid 


halo 


89 


0.01 85 


10 








1 




spheroid 


8 


2.5 40 


1.5 


14 






0.6 




disk 


11 


0.6 40 




18 


7 


0.3 




Massive disk 


halo 


79 


0.01 110 


8 








1 




spheroid 


4.4 


4.5 40 


1 


9 






0.6 




disk 


19 


1.4 40 




26 


6 


0.3 




Reference 


halo 


123 


0.065 100 


4 








1 




spheroid 


4.4 


4.5 40 


1 


9 






0.6 




disk 


5.3 


0.35 40 




8 


6.4 


0.3 





We do not include any nuclear bulge component mainly because 
the bulge region is unlikely to be useful to experiments wishing to 
measure any near-far spatial asymmetry arising from Machos, but 
also because not including it allows us to test the widest possible 
mass range for the other components. 

Constraints on possible models come from the M31 surface 
brightness profile, the rotation curve and from dynamical mass es- 
timates, so we ensure that our models satisfy at least these con- 
straints. We also ensure that the B-band M/L values of the visible 
populations in our models is not lower than considered reasonable 
for those populations, though we do not impose any upper bound 
on M/ Lb- 

For our halo we assume a simple cored near-isothermal 
sphere, with a density falling off with radius r as 

Po 



Phalo 



l + (r/a)2' 



(5) 



where po is the central density and a is the core radius. The disk is 
assumed to have a sech-squared distribution: 



Pdisk = po exp sech^ ^-^^ 



(6) 



where (a, z) are cylindrical coordinates, h the disk scale length and 
H the scale height. Finally, we model the spheroid by a flattened 
power law: 

Po 



Psph = 



1 + [(cr/a)2 +g-2(z/a)2]»/2' 



(7) 



with q the ratio of minor to major axes an d a the power law \n- 
dex. iReitzel. Guhathakurta & Gould fl99? ) find that s ~ 3.8 and 
g ~ 0.6 provide a good fit to their number counts of RGB stars at a 
projected distance R = 19 kpc from the M31 cent r e, toge ther with 
number counts from iHoUand. Fahlman & Riched il996h in fields 
at J? = 7.6 and 10.8 kpc. However, a single power law does not 
provide a good global description of th e M31 spheroid, which is 
better fit by a de Vaucouleurs profile iPritchet & van den Berghl 
Il994h . However, as we shall see, this simple parameterisation re- 
produces the surface brightness distribution remarkably well in the 
inner galaxy. The main effect for microlensing calculations in us- 
ing a single power law rather than a de Vaucouleurs profile is to 
increase the spatial dispersion of spheroid stars along the line of 
sight and so increase their contribution to the microlensing rate. 
For such a steep profile this is a small effect and, in any case, pro- 



vides for a worst-case assessment of the impact of spheroid stars 
for Macho microlensing surveys. 

The parameters of the galactic components are listed in Ta- 
ble Q for the three extreme models. The massive halo model as- 
sumes an extreme total M31 mass of around 2 x 10^^ M!0> consis- 
tent with traditional estimates, whilst the other two models assume 
a total mass closer to 1 x 10^^ M^t^, in line with recent determina- 
tions iEvans & WilkinsonI EoOOt iGuhathakurta. Reitzel & Grebell 
l200Cl) . The M/Lb values for the disk and spheroid in the mas- 
sive halo model are close to minimal and do not allow much 
room for non-luminous matter. By contrast in the massive spheroid 
model the disk and spheroid have larger M/Lb values. This 
model possesses a halo with a large core radius and a disk with 
a slightly larger scale length in order to increase the spheroid 
contribution. In the massive disk and spheroid models the high 
disk M/Lb assumes a dark matter dominated disk, something 
which is known not to be the case for the Milky Way disk 
iCreze. Chereul. Bienavme & Pichonll998l) . The disk and spheroid 
M/Lb values are comparable for the massive spheroid model, 
whilst the spheroid M/Lb is higher than that of the disk for the 
massive halo model and lower for the massive disk model. We have 
also included a reference model in TableQ This model employs a 
less extreme set of parameters and therefore can be considered a 
more likely candidate model than the extreme models. 

The resulting surface brightness profiles and rotation curves 
of the models are plotted in Figure □ The left-hand fi gures show 
the predicted surface brightness distribution for the models, in- 
cluding the individual contributions of the spheroid and disk. Also 
plotted on these figures are the radially-averaged measurements of 
IWalterbos & Kennicutll ^9^). All models provide a remarkably 
good fit to the data. The right-hand plots show the rotation curve 
contributions for the spheroid, disk and halo, along with the over- 
all profile. Rotation curve data, taken from Kent ( 1989), are also 
displayed. Whilst there are differences in the rotation curves of the 
models, they provide reasonable fits to the data given their simplic- 
ity. 

We must emphasize that the assumed spheroid density law 
in equation Q is unlikely to provide an accurate description to- 
wards the inner galaxy, despite the good agreement with the surface 
brightness profile in Figure Q. However, this is not particularly a 
problem for our stu dy of the distribution of microlensing events 
in the outer regions. iReitzel. Guhathakurta & Gould jl998l> found 
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Figure 1. The _B-band M3 1 major axis surface brightness profile (left) and rotation curve (right) for the massive halo model, massive spheroid model, massive 
disk model, and reference model. The individual contributions of the spheroid, disk and halo are shown by the light dot-dashed lin e, medium d ashed line and 
dark dotted line, respectiv ely, whilst th e overall profile is shown by the black solid line. The surface brightness measurements of lWalterbos & Kennicuti 1.19871) 
and rotation curve data o f lKenlHl989t) are shown by the crosses. 
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that the spheroid luminosity density in M3 1 is about ten times that 
of the Milky Way at a location equivalent to the Solar position in 
the Milky Way. For the spheroid model in the massive halo and disk 
models, which employs a standard spheroid M/ Lb value, we set 
the mass density normalisation to be consistent with this ratio. For 
the massive spheroid model the mass density is twice as large due 
to the larger assumed M/ L b ■ 

It should be stressed that the assumed form of the populations, 
especially that of the halo, plays a crucial role, along with the ro- 
tation curve and surface brightness distribution, in determining the 
relative masses of the components within each model. For exam- 
ple, in the massive halo model there is some leeway in the allowed 
ratio of spheroid to disk mass which would still produce acceptable 
profiles. Where this is the case we have tended to maximise the 
spheroid contribution. However, if we were to allow for dark halo 
distributions with density profiles which could be steeper or shal- 
lower than r^^ then obviously an even greater range of spheroid 
and disk mass ratios would be permissible. Our conclusions must 
therefore be viewed within the context of our model assumptions. 



4 MICROLENSING EVENT MAPS 

Whilst our three-component models satisfy the observed rotation 
curve and surface brightness pro file, we must remain aware of 
the limitations of their simplicitv. IWidrow. Perrett'& 
have undertaken a detailed investigation of a number of self- 
consistent three-component bulge-disk-halo models of M31 em- 
ploying a range of additional constraints including the inner ve- 
locity dispersion, bulge and disk MjL, dynamical mass estimates, 
disk stability and cosmological considerations. Such analyses are 
especially important in constraining the expected rate of microlens- 
ing. However, we can probe the relative contributions of each pop- 
ulation within our models by instead considering the optical depth, 
r, a quantity which depends only upon the density distribution of 
lenses and sources. In this sense it is a more robust quantity than 
the event rate which depends additionally upon the lens and source 
velocity distributions as well as the lens mass function. Physically, 
T is the number of ongoing microlensing events per source star at 
any instant in time. So, we can easily compute the instantaneous 
event number density as a function of position by multiplying r by 
the number density of sources. 



4.1 Pixel lensing event density and asymmetry 

Pixel lensing differs from conventional microlensing in that the 
source stars are mostly unresolved, so microlensing events are seen 
only as small enhancements of the local surface brightness flux. 
The spatial distribution of pixel lensing events therefore depends 
upon the source luminosity function and the background surface 
brightness distribution, in addition to the usual parameters of clas- 
sical microlensing. Specifically, a sufficiently high event magnifi- 
cation requires that the impact parameter, u, in units of the Einstein 
radius be less than some threshold, ut, where ut depends upon 
the unlensed source flux F and the background surface brightness 
at the event location. The number of microlensing events which are 
detectable at some sky position (x, y) will therefore depend upon 
the quantity 



where (f> is the source luminosity function. Since the optical depth 
of events with impact parameter less than u scales as , we define 
the pixel-lensing optical depth to be 



_2 



(9) 



UT{x,y) 



J UT{F,x,y)4>{F)dF 



(8) 



This quantity is the number of ongoing events with u < ut at a 
given location. 

We compute Tp over a grid of sky positions {x,y) on the M31 
disk for each of the models in TableQ and for each combination of 
lens and source population. The luminosity function of M3 1 stars 
is known only for the brightest stars so, to calculate ut, we make 
l/-band calculations taking M31 sources to have a Solar neigh- 
bourh ood luminosity function for My > lBahc air& Soneir3 
ll98CtlWielen. Jahreiss & KrugeJl987h and use the M31 luminos- 
ity function determined bv lHoUand^ahlm an & Rich^ il996h for 
Mv ^ 0. Whether a microlensing-induced flux enhancement is de- 
tectable depends upon experimental sensitivity, so Tp is implicitly 
an experiment-specific quantity. We base our compu tations loosely 
on the POINT-AGAPE survey l lAuriere et ai]|200ll) by computing 
flux changes within a square pixel array (a "super-pixel") with a 
2'.'3 field of view. The radially-averaged surface brightness mea- 
surements of W alterbos & Kennicutt (,1987.) are used to compute 
the M31 surface brightness within the super-pixel. We assume 
that typically 40% of a randomly-positioned point source PSF 
is contained within the super-pixel and that the survey is sensi- 
tive to super-pixel flux variations exceeding 1%. In addition to 
the M31 surface brightness, we assume a sky background in V 
of 19.5 mag/arcsec'^, which conservatively allows for a lunar sky 
background contribution. 

Table|2|displays the total number of ongoing events over a cen- 
tral 100 X 70 arcmin^ region of M31, within which current surveys 
are looking for Machos. We have assumed for this Table that the 
M3 1 halo is full of Machos; for a partial Macho halo the number 
of halo events scales linearly with halo fraction. Since we are in- 
terested in the effect of stellar lenses on the spatially extended Ma- 
cho population we have defined an 8 arcmin radius exclusion zone, 
counting only events outside this region. Since the inner region 
may be dominated by stellar lensing, surveys are likely to define 
such masks as a simple way of minimising the impact of the stellar 
lens contribution, so making their results less sensitive to the pre- 
cise details of the inner galactic structure. The left-hand columns 
of the Table shows the number of ongoing events for each popu- 
lation, together with the cumulative total. The massive halo model 
provides the largest number of events overall, thanks to the sub- 
stantially more massive halo assumed for this model. For all three 
extreme models M3I Machos provide the largest contribution of 
any individual lens population. However, whilst M3 1 Machos con- 
tribute almost 90% of the event density for the massive halo model, 
they provide only around half of the events for the massive spheroid 
and disk models. In the far disk, beyond the region considered, the 
relative contribution of M31 Machos is significantly larger but the 
number of available sources is few and these regions are not being 
targeted by current surveys. 

The disk shows a factor of 6 variation in the expected num- 
ber of events between the three extreme models, whilst for the halo 
and spheroid the variation is just over a factor of 2. The variation 
in the overall number of events is rather less, at around 40%. As 
we shall discuss in Section 1421 the number of events computed 
from Tp does not necessarily reflect the relative contributions one 
might observe in the event rate, since the rate also depends upon 
the event durations which are sensitive to the lens and source ve- 
locity distributions, as well as the lens mass function. Differences 
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Table 2. The number of ongoing events and the number asymmetry for the three extreme M31 models and for the reference model. The 
asymmetry is expressed both in absolute numbers of events and as a ratio (in brackets). The figures assume a full Macho halo and apply to 
the central 100 X 70 arcmin^ region orientated along the major axis of M31 (see Figure|2|. We exclude events occurring within a radius of 
8 arcmin of the M3 1 centre. 



Model 


No. of events (100 X 70 
Halo Spheroid Disk 


arcmin^ ) 
Overall 


Stellar 


Asymmetry (far 
Halo 


disk/near disk) 
Halo-to-disk 


Overall 


Massive halo 


1.27 


0.12 


0.07 


1.46 


0.12/0.07 (1.67) 


1.01/0.26 (3.92) 


0.83/0.16(5.24) 


1.13/0.33 (3.43) 


Massive spheroid 


0.70 


0.30 


0.27 


1.27 


0.29/0.28 (1.02) 


0.51/0.19 (2.75) 


0.37/0.09 (3.87) 


0.80/0.47(1.72) 


Massive disk 


0.51 


0.12 


0.41 


1.04 


0.21/0.32 (0.68) 


0.38/0.13 (3.07) 


0.32/0.08 (3.77) 


0.60/0.44(1.36) 


Reference 


1.08 


0.11 


0.11 


1.31 


0.12/0.10(1.22) 


0.85/0.23 (3.60) 


0.69/0.14(4.80) 


0.97/0.34 (2.89) 



in these distributions between the lens populations will therefore 
mean differences between the relative contribution to the rate and 
the contribution to Tp. 

The right-hand columns of Table|2|show the asymmetry in the 
number of ongoing events against the far disk to those against the 
near disk. The numbers for the halo population again assume a halo 
full of Machos. The signature of asymmetry is the key diagnostic 
for the identification of Machos, but as Table |2| shows, one may 
also see an asymmetry in spatial distribution of stellar events. In 
the massive halo model the ratio of far-disk to near-disk events is 
3.9 for M31 Machos and 1.7 for stellar lenses. Stellar lenses in this 
model therefore have a spatial asymmetry which is weaker but in 
the same direction as that of M31 Machos. For the massive disk 
model the stellar asymmetry is reversed, with only 40% of stel- 
lar events occurring on the far side of the disk. This is because 
the disk in this model has a very high M/ L and so disk lensing 
against spheroid sources dominates the number of stellar events. 
In the massive spheroid model the M/L of the spheroid and disk 
components are comparable and the distribution of stellar events is 
close to symmetric. These results bear out the analytic prediction 
of Section |2| that the presence and direction of spatial asymmetry 
due to stellar lenses is dictated by the ratio of M/L values in the 
spheroid and disk populations. Clearly, unless there is fine tuning in 
the spheroid and disk M/ L, one should expect some asymmetry in 
the spatial distribution of stellar lenses. However, this asymmetry 
should be weaker than that due to M3 1 Machos if Machos con- 
tribute significantly to the dark matter. 

The column labelled "Halo-to-disk" in Table |2| shows the 
asymmetry of M31 Machos against disk sources alone, ignoring the 
effect of spheroid sources. Comparison with the adjacent column 
labelled "Halo", which includes both disk and spheroid sources, 
indicates that the presence of the spheroid sources dilutes the M3I 
Macho asymmetry by 20 — 30% for the three extreme models. This 
result is largely insensitive to the very different mass normalisa- 
tions assumed for the three extreme models. 

In Figure|2|the left-hand panels display contours of event num- 
ber density, that is Tp multiplied by source number density, for both 
M3 1 Macho and stellar lensing events over a central area covering 
100 X 70 arcmin^. A full Macho halo is assumed. The x and y 
axes are orientated along the major and minor axes of M31, respec- 
tively, with positive y pointing towards the near side of the disk. 
The dark contours show the event density of M3 1 Machos, with the 
thick contour indicating a density of 3 x 10~^ events/arcmin^. The 
density increases towards the M3 1 centre by a factor of 2 for each 
successive contour. The near-far asymmetry signature of the M3I 
Macho population is strongest for the massive halo model (top) but 
is clear for all three extreme models. The different mass normalisa- 
tions of each component in the three extreme models have a clear 



impact on the number of M31 Machos, as evidenced by the relative 
positions of the thick contours. 

The stellar lens distribution is shown by the light grey con- 
tours in the left-hand panels of Figure |2| Again, the thick con- 
tour indicates an event density of 3 x 10^* events/arcmin^ and 
contour spacing is the same as for the Macho distribution. As ex- 
pected, the contours for the massive disk model are substantially 
flattened compared to those for the massive spheroid model. The 
stellar asymmetry for the massive halo and massive disk models 
is clearly evident. In the massive halo model the M3 1 Macho dis- 
tribution is clearly more extended than that of the stellar lenses, 
however for the massive spheroid and disk models the stellar con- 
tribution remains substantial out to large distances from the M3I 
major axis. 

The relative importance of the stellar and Macho contribu- 
tions to the number of ongoing events varies considerably across 
the three extreme models, as indicated by the right-hand panels 
of Figure |2| The light grey contours show the region where stel- 
lar lensing dominates, whilst darker contours indicate areas where 
M31 Machos provide the dominant contribution, assuming a full 
Macho halo. The locus of positions where the stellar and M3 1 Ma- 
cho event number densities are equal is shown by the dotted line. 
Inside the dotted line the ratio of stellar-to-Macho event densities 
doubles as one progresses towards the M31 centre, whilst outside 
the dotted line the stellar-to-Macho event ratio halves for each con- 
tour level going outwards. The positions of the INT fields surveyed 
by the POI NT-AGAPE and MEGA teams are al so indicated for 
comparison jAuriere et alj200llJCrotts et alj200lh . In the inner re- 
gions the event density is dominated by M31 Machos only for the 
massive halo model (top). In the massive spheroid model (middle 
panel) stellar lenses dominate within the central 10 — 15 arcmins 
and also over a large area of the near disk between | a; | < 30 arcmins 
and < y < 20 arcmins. For the massive disk model stellar events 
provide the bulk of the near-disk events over almost the whole re- 
gion considered. However, even for these extreme models Machos 
dominate the far disk event numbers for the case of a full Macho 
halo. For a halo with only a 50% halo fraction the dotted contour is 
shifted outward by one contour level and for a 25% halo fraction it 
is shifted outward by two levels. 

From this we can estimate how small the M31 Macho halo 
fraction would have to be for the POINT-AGAPE/MEGA INT 
fields to be dominated everywhere by stellar rather than Macho 
lensing events. For the massive halo model Machos should dom- 
inate at least part of these survey fields provided their halo fraction 
exceeds 2%. For the massive disk model we require a halo fraction 
above 4% and for the massive spheroid model the Macho fraction 
must be at least 8%. More generally, at least half of the far-disk 
region considered must be dominated by Machos if the M3 1 Ma- 
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Figure 2. M31 contour maps of instantaneous pixel-lensing event number density (left) and of relative Macho-to-steUar event nvmiber density (right) for 
the massive halo, massive spheroid, massive disk and reference models, assuming a full Macho halo. The total event nvmiber density contribution due to 
M31 Macho lenses is shown by the dark contours in the left-hand panels, whilst the lighter contours show the contribution due to stellar lenses. The thick 

contours represent a density of 3 X 10^'' events/arcmin^ . Successive contours represent a factor 2 change in event density. The ratio of the Macho and stellar 
contributions is shown in the right-hand panels. The dark contours show regions where M31 Macho lenses dominate, whilst lighter contours show where 
stellar lensing dominates. The black dotted contour shows where the M31 Macho and stellar lens contributions are equal. Again, successive contours represent 
a factor 2 ratio change. The tield positions of the POINT-AGAPE/MEGA INT survey are also shown. 
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cho halo fraction exceeds about 25%. Interestingly, though stellar 
lenses dominate the largest fraction of the near disk in the massive 
disk model, their density drops off faster in the far disk than in the 
massive spheroid model, which is why the spheroid model requires 
a larger Macho halo fraction. 

The event density distribution for the reference model is also 
shown in Figure |2| The distribution for this model is intermediate 
to those of the massive halo and spheroid models. The event den- 
sity is dominated by Machos for this model outside of the innner 
5 arcmins, assuming a full Macho halo, and dominates both the 
near and far disk. In this respect, the massive disk model may be 
regarded as the most extreme and unlikely of our three extreme 
models. 



4.2 Event rates and timescales in pixel lensing 

The pixel-lensing optical depth, Tp, can be related to pixel-lensing 
timescales and rates in much the same way as the analogous ob- 



servables in classical microlensing 


;. We can define the pixel-lensing 


rate as Fp = utV i 


Kerins et al. 


2001), where F is the classical 



microlensing rate. Whilst in classical microlensing the character- 
istic event duration is given by the Einstein Radius crossing time, 
tE, in pixel lensing iE is generally not measurable because of a 
near-degeneracy in the light-curve between the event timescale, 
magnification and source luminosity (Gould .1996.) . Instead the 
usual timescale observable in pixel lensing is the full-width half- 
maximum duration of the event. This can attain a maximum given 
by 



2^U'TtE T <. 1) 

ITdaysfi^H^) 



''1/2 



, 100 km s 

where m is the lens mass, v is the lens transverse velocity relative 
to the line of sight and D ~ Ds — Di. With these definitions we 
have 



TT _ 

T-o — ^tl/2Fp, 



(11) 



where Ii/2 — 2\/Z%li^ is the average ti/2 timescale and over-lines 
on other quantities also denote averages. Equation <1 H relies on the 
fact that u = ut/2, so ti/2 — t^f^ 1"^ average for fixed tE. 

We can estimate the typical t^j^. duration of pixel-lensing 
events from very general considerations of the various lens and 
source components. Though details of the lens and source velocity 
distributions can vary even for distribution functions with similar 
density profiles, the event duration represents a line-of-sight statis- 
tical average, and so is a reasonably robust quantity. Furthermore, 
the duration depends not just on the lens and source distribu- 
tion functions but also on the lens mass function and source lumi- 
nosity function. The timescale distribution is therefore rather 
broad and variations in the lens or source velocity distribution make 
a relatively modest contribution to this width. 

The underlying tE depends upon the lens mass m, relative 
transverse velocity v and lens-source separation D. Additionally, 
ti/2 depends upon the impact parameter u < mt, where ut de- 
pends upon the source luminosity and background surface bright- 
ness. For our computations of Tp in the previous subsection, ut ~ 
0.1 over most of the area considered. On average we should find 
u = ut /2 ^ 0.05, so we adopt this value for all estimates. 



The typical mass of stellar lenses depends on their mass function 
(f>{m). In the Solar neighbourh ood cf) oc m"^'^^ for m > 0. 5 Mq 
and m^i * below 0.5 Mq IGould. BahcaU & Flvnrj[l997h . The 
lensing rate per unit lens mass is dV /dm oc m}^^<j>{m), and there- 
fore the mean lens mass for stars between 0.01 Mq and 10 Mq 
is m ~ 0.5 Mq in the Solar neighbourhood. In the absence of 
detailed knowledge of the M31 stellar mass function we adopt 
m — 0.5 Mq for M31 stars. This also corresponds to the mass 
scale for Machos favour ed by the resul ts of the MACHO microlens- 
ing survey of the LMC jAlcocketall2 000). 

The relative transverse velocity v varies according to lens and 
source population. If lens and source velocities can be described 
by a combination of rotational and random motions then we can 
expect the characteristic velocity to be u ~ [(t'sys.i — t^sys,s)^ + 
'''rani + '»Ta.ii,sY^^, whcrc subscopt "1" and "s" denote lens and 
source velocities, respectively, "sys" denotes systematic motions 
due to rotation and "ran" represents the random velocity compo- 
nent. Table|3|lists characteristic values for v for each combination 
of lens and source population. The numbers assume that Machos 
have 2-D random velocities of 235 km s^^, comparable to the disk 
rotation speed. We also assume disk and spheroid 2-D random ve- 
locities of 45 and 170 km s^^, respectively, and a spheroid rota- 
tional component of 30 km s~i. The contribution of rotational mo- 
tions to V will depend upon where events occur against the inclined 
M31 disk but, on average over the disk, the transverse component 
is about 80% of the total rotational motion. Table|3]also indicates 
typical values for the lens-source separation D. 

Estimates for the typical ti/2 and Fp are computed from these 
representative values of m, v and D, together with an assumed im- 
pact parameter u — 0.05. Based on these numbers, for disk lenses 
we can expect TV ~ 8 events/yr for the massive halo model, or up 
to six times this number in the case of the massive disk model. For 
spheroid lenses we can expect ~ 18 events/yr from the massive 
halo and disk models ranging up to 43 events/yr for the massive 
spheroid model. The expected number of M31 Machos ranges from 
~ 32 events/yr for the massive disk model up to ~ 80 events/yr for 
the massive halo model. 

The ratio of halo:spheroid:disk event rates is 1:0.2:0.1 for the 
massive halo model, 1:1:0.75 for the massive spheroid model and 
1:0.6:1.4 for the massive disk model. These numbers indicate that 
in the massive disk and spheroid models stellar events are expected 
to dominate the overall rate. This is because the disk and spheroid 
in these models are effectively repositories of large amounts of dark 
matter in their own right. For the massive halo model, which em- 
ploys more conventional values for the disk and spheroid M / L, 
Machos out-number stellar events by more than 3:1, with Machos 
becoming increasingly dominant in the far disk. 

An important caveat to these estimates is that they assume a 
100% detection efficiency for events of all timescales, something 
which is not the case for any real experiment. Sampling efficiency 
in particular is likely to have a large impact on the detected num- 
ber of stellar events, which for the most part have typical durations 
of just a few days. Current experiments are unlikely to have good 
detection efficiency for events with durations below ti/2 ^ 3 days. 
Given that the log(ti/2) distribution is expected to be very broad, 
and the fact that we estimate Ii/2 ~ 3 days for most stellar events, 
we can expect that many more stellar events will escape detection, 
as compared to Machos with m ~ 0.5 Mq. This may mean that, 
even for the massive disk and spheroid models above. Machos still 
dominate the rate of detected events. Additionally, finite source size 
effects are not taken into account in these estimates, but are likely 
to be more important for stellar lenses than for Machos, due to 
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Table 3. The estimated event rate, assuming 100% detection efficiency, for tlie tliree extreme 
M31 models. The figures assume a full Macho halo and apply to the central 100 X 70 arcmin-^ 
region of M3 1 , though excluding events within a radius of 8 arcmin of the M3 1 centre . D, 
V and ti^2 indicative values, N{Tp) is the number of ongoing events computed from Tp 
and N(rp) is the estimated event rate, taking the indicative ti/2 t" be an estimate of 
A lens mass of 0.5 Mq is assumed for ti/2 ^nd Fp. The numbers in brackets correspond to 
predictions for the massive (halo, spheroid, disk) models, respectively. 



Lens/Source 


D/kpc 


ti/km s ^ 


ti/2/days 




Af(rp)/year 


halo/disk 


10 


300 


6 


(1.24,0.56,0.50) 


(63,29,25) 


halo/spheroid 


10 


290 


7 


(0.28,0.25,0.11) 


(17,15,7) 


spheroid/disk 


1 


240 


2 


(0.12,0.25,0.12) 


(15,29,15) 


spheroid/spheroid 


2 


240 


4 


(0.03,0.12,0.03) 


(3,14,3) 


disk/disk 


1 


60 


10 


(0.03,0.08,0.19) 


(1,3,8) 


disk/spheroid 


1 


240 


2 


(0.05,0.26,0.30) 


(7,30,36) 


TOTAL 








(1.75,1.52,1.25) 


(106,121,94) 



the smaller lens-source separation. When D is sufficiently small 
that the Einstein radius has a comparable angular size to the source 
star one expects differential magnification across the face of the 
source which may make stellar lens events harder to detect. Rela- 
tive rates may also be affected by the orientation of Macho orbits. 
The Macho velocities in Table|3|assume an isotropic velocity distri- 
bution, whereas Machos on radial orbits are expected to produce a 
larger event rate, and an enhanced near-far asymmetry iKerins et alJ 

Perhaps surprisingly, the total number of events predicted for 
the three models is similar. The ratio of the total event rates for the 
massive halo, spheroid and disk models is 1:1.1:0.9, respectively, 
slightly different from the 1:0.9:0.7 ratio for the total numbers of 
ongoing events computed from Tp. The difference is due mainly to 
the fact that stellar lenses typically have smaller tip durations, so 
models with larger stellar lens contributions show an increase in 
the rate relative to their optical depth contribution. Overall the total 
number of events is a remarkably robust quantity given the large 
differences between the models. 



5 DISCUSSION 

One of the key differences between M31 and our own Galaxy is 
that M31 possesses a more massive and luminous stellar spheroid. 
We have assessed the consequence of this for pixel-lensing exper- 
iments which are targeting M3 1 . These experiments are hoping to 
confirm the existence of massive compact halo objects (Machos) 
by detecting an excess of events towards the far side of the M31 
disk. The near-far spatial asymmetry expected for a spheroidal dis- 
tribution of M31 Machos is regarded as a crucial test of the Macho 
hypothesis. 

In this study we have shown that the presence of a massive 
stellar spheroid can induce asymmetry in the spatial distribution of 
stellar lenses. Analytical predictions, confirmed by more detailed 
model calculations, show that stellar lenses in the disk and spheroid 
are asymmetrically distributed if the spheroid and disk mass-to- 
light (M/L) ratios differ. A higher spheroid AI/L should result in 
a stellar lens asymmetry in the same direction as that of Machos 
(i.e. an excess of far-disk events), whilst a higher disk M / L re- 
sults in an asymmetry of opposite sign to the Macho asymmetry 
(i.e. an excess of near-disk events). A positive consequence of this 
is that two competing dark matter scenarios, one in which one has 
a spheroidal dark Macho halo and one in which the dark lenses oc- 



cupy a massive disk or highly flattened halo, should be distinguish- 
able by the fact that they predict spatial asymmetries of opposite 
sign. 

We have constructed three extreme models of M3 1, each com- 
prising disk, spheroid and Macho halo populations. The massive 
halo model assumes conventional AI/L values for the spheroid 
and disk whilst the massive spheroid and disk models assume high 
AI/L values consistent with dark matter dominated stellar popu- 
lations. By defining the pixel-lensing optical depth we have com- 
puted the number of ongoing events for each population and have 
constructed microlensing maps over a region of M3 1 currently be- 
ing monitored by a number of survey teams. Whilst Machos make 
up 90% of the events in the massive halo model, stellar lenses pro- 
vide at least half of the signal for the other two models, even for a 
full Macho halo. However, even for these extreme stellar models, 
M3 1 Machos should dominate over most of the far disk if their halo 
fraction is around 25% or larger. 

We find for all three models that the inclusion of the spheroid 
dilutes the Macho near-far asymmetry signature by about 20 — 30% 
over the disk-only asymmetry values. The overall asymmetry varies 
from 1.4 and 1.7 for the massive disk and spheroid models to 3.4 
for the massive halo model. Kerins et al. (2003) argue that asym- 
metries weaker than a factor of 2 will be difficult to detect with 
current surveys, so if the M31 disk and spheroid is strongly dark 
matter dominated (in contrast to our own Galaxy) they could con- 
ceivably mask the asymmetry signature due to Machos. However, 
from general considerations of the likely Macho and stellar lens ve- 
locities we argue that many stellar lenses have durations too short to 
be detected by the survey teams. If Machos have a mass of around 
0.5 M ^, as favoured by the MACHO LMC survey I Alcock et a^ 
'200(f) then the overall asymmetry of detected events is likely to be 
significantly larger than these estimates because the longer dura- 
tion Macho events are easier to detect. Estimates of the expected 
event rates indicate that the total number of pixel-lensing events 
is remarkably insensitive to the relative mass normalisation of the 
halo, spheroid and disk components, though, again, there is likely 
to be a strong bias towards the detection of stellar-mass Machos 
over stellar lenses because of the short duration of stellar lensing 
events. 

In summary, the presence of the spheroid complicates matters 
for pixel-lensing surveys of M31 and cannot be ignored. A pos- 
itive detection of near-far asymmetry cannot, by itself, be taken 
as a confirmation of the Macho scenario since a stellar spheroid 
will produce a similar result if its AI/L is larger than that of the 
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disk. In this case an analysis of the spatial extent is likely to be 
able to distinguish Machos from spheroid lenses, as spheroid lenses 
should be more strongly concentrated toward the M31 centre. The 
spheroid also has beneficial consequences for the microlensing sur- 
vey teams. Candidate samples may be dominated by Machos, or by 
lenses in a highly flattened halo or disk, or they may instead be ho- 
plessly contaminated by variable stars mistaken for microlensing. 
In the absence of a spheroid it might prove very difficult to distin- 
guish between the latter two cases. The presence of the spheroid 
provides an all-important third signature with which one may be 
able to distinguish between all three scenarios. In the Macho case 
one expects a near-far asymmetry. In the case of a highly flattened 
lens population one expects a far-near asymmetry (i.e. an asymme- 
try of opposite sign to that of Machos) due to the presence of the 
spheroid. Finally, if the candidates are mostly variable stars their 
distribution should be S3fmmetric. 
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